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The spreading of an incompressible viscous liquid over an isotropic homogeneous unsaturated porous substrate 
is considered. It is shown that, unlike the dynamic wetting of an impermeable solid substrate, where the 
dynamic contact angle has to be specified as a boundary condition in terms of the wetting velocity and other 
flow characteristics, the 'effective' dynamic contact angle on an unsaturated porous substrate is completely 
determined by the requirement of existence of a solution, i.e. the absence of a nonintegrable singularity in 
the spreading fluid's pressure at the 'effective' contact line. The obtained velocity dependence of the 'effective' 
contact angle determines the critical point at which a transition to a different flow regime takes place, where 
the fluid above the substrate stops spreading whereas the wetting front inside it continues to propagate. 



1. Introduction 



The modelling of the spreading of liquids over porous substrates in the framework of continuum mechanics 
requires, and is based on, th e separation of scales be tween the 'macroscopic' (or 'Darcy-scale') and 'microscopic' 



(or 'pore-scale') processes ([Barenblatt et 



1990). In the continuum approximation, the description of the 
spreading phenomenon brings in the notion of an 'effective' smooth penetrable solid substrate, which is how 
the actual porous medium is represented, together with the notions of 'effective' contact lines and 'effective' 
contact angles that the free surface of the pure fluid above and the wetting front insid e the substrate form 



1999), are 



with it. These notions, being averages in the sense of mechanics of multiphase media ([Whitakerl 
fundamentally different from the concepts of 'contact line' and 'contact angle' used in the modelling of dynamic 
wetting on the actual, as opposed to 'effective', impermeable solid surfaces. In particular, the notion of a 'static' 

(or 'equilibrium') contact angle, central to the modelling of dynamic wetting where it is a measure of wettability 
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of the solid (Ralston et al 
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20081 ). becomes mean ingless for the Darcy-scale description of the liquid spreadin g 



over a porous medium: as experiments show (e.g. 



Clarke et al. 


2002 


Starov et al. 


2003 


Markicevic et al. 


2010) 



if, say, a drop of liquid is deposited onto an unsaturated wettable porous substrate, the eventual equilibrium 
state will be the drop completely imbibed into the porous medium, with no liquid left above it and hence no 
'effective' static (equilibrium) contact angle between the, now non-existent, free surface of the pure liquid and 
the effective surface of the substrate. 

The main theoretical implication of this absence of a meaningful effective static (equilibrium) contact angle 
is that the effective dynamic contact angle that the free surface of a spreading liquid forms with a porous 
substrate on the Darcy-scale cannot be regarded as essentially a perturbation of the static contact angle, which 



Dussan V. 


1979 


de Gennes 


1985 


Blake 


2006 


Shikhmurzaev 


2011 



for reviews). Therefore, it becomes important 



to consider the problem from first principles, without implying a-priori that the case of a porous substrate can 
be described by adjusting the concepts borrowed from the modelling of dynamic wetting of impermeable solid 
surfaces, such as the equilibrium contact angle. 

In the present paper, we show that, unlike the situation one has in dynamic wetting, where the 'microscopic' 
dynamic contact angle has to be specified as an additional boundary conditioij~j, for the process of the fluid 
spreading over a porous substrate one can find the 'effective' dynamic contact angle from the requirement of 
the absence of a nonintegrable singularity of the fluid's pressure at the contact line, i.e. essentially from the 
requirement that a solution exists. (The flow inside the porous matrix is treated in the standard way, i.e. 
employing a dynamic wetting model for the impermeable solid surface to describe the propagation of menisci 
in the pores.) The limits confining the considered regime suggest certain experimentally verifiable predictions 
of the model which we briefly discuss in the light of the available experimental data. 



2. Problem formulation 

Consider an incompressible Newtonian fluid of density p and viscosity /i spreading at a speed U over an 

isotropic homogeneous unsaturated porous substrate characterized by an effective pore size a. The gas displaced 

by the fluid from the substrate and from the inside of the porous matrix is assumed to be ideal and at a constant 
f Here we are not discussing the so-called 'apparent' contact angle resulting from the free-surface bending near the 

contact line. This angle is not part of the mathematical problem formulation; it is an auxiliary concept introduced in 

some works and in different ways (see, e.g., ?) to interpret experimental results when the measured angle differs from 

the angle imposed in the problem formulation. 
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Figure 1. Sketch of the spreading of a viscous liquid over an unsaturated wettable porous substrate in a coordinate 
frame moving with the contact line in the framework of continuum mechanics. The magnified view illustrates the main 
flow mode on the pore scale, where 8d is the dynamic contact angle formed by the moving meniscus with the wall of the 
pore. Regions fii, Q.2 and correspond to the pure liquid, the saturated part of the porous medium and the unsaturated 
porous matrix, respectively; 9d is the 'effective' dynamic contact angle formed by the free surface of the liquid Sf with 
the surface So of the 'effective' porous substrate, and S w is the wetting front inside the substrate. 

pressure with respect to which the pressure in the pure fluid and in the fluid inside the porous medium will 
be measured. In order to be able to model the process in the framework of continuum mechanics, we need 
a separation of scales between the macroscopic and the pore-scale phenomena, i.e. we have to consider the 
continuum limit 



where L is the characteristic length scale on which the phenomenon is described. The resulting model will 
be applicable to experiments if a/L <C 1, with the actual value of a/L determining its accuracy. In the Oth 
approximation in the above limit, one has (see Fig.[I| (a) a macroscopic 'wetting front' as a sharp interface S w 
separating the saturated porous medium from the unsaturated matrix ^3, (b) an effective 'contact line' at 
which the free surface Sf, confining the domain Qi occupied by the pure fluid, and the wetting front meet, and 
(c) well-defined 'contact angles' 8rj and 9 W that the free surface and the wetting front form with the 'effective' 
surface So of the solid. The reference frame in which the problem will be considered and the directions of unit 
normals n, n and no to, respectively, S w , Sf and Sq are shown in the figure. Importantly, in the scheme outlined 
above we have already made an assumption that the two contact lines, i.e. the contact line CL1 formed with the 
substrate by the free surface Sf and the contact line CL2 formed with it by the wetting front S w , coincide. This 
is always what happens when the fluid is first brought in contact with the substrate, and we will examine what 
follows from this initial situation and, later, the conditions when the assumption that CL1 and CL2 coincide 
no longer holds. 

We will consider the case of small Reynolds and capillary numbers for the flow of the pure fluid, i.e. the limit 




(2.1) 



Re = pLU I [i — ► and 
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a 



(2.2) 



where a is the surface tension of the fluid-gas interface. In this limit, to leading order, inertial effects can be 



neglected, and, as in the case of the dynamic wetting of an impermeable solid surface (|Huh fc Scriven 



S hikhmurzaev 



1971 



19931 ). from the normal-stress boundary condition on the free surface we have that near the 
contact line the free surface (Sf) is locally planar, so that locally Oi is a wedge region. For simplicity we will 
also neglect gravity, though its inclusion would not change the main results of the analysis below. 

Scaling the pressure and velocity in fli with fiU / L and U respectively, one has that in Ox the dimensionless 
pressure p and velocity u obey the Stokes equations 



V-u = 0, Vp = V 2 u, (refit) 



and on the free surface satisfy the standard kinematic and tangential-stress boundary conditions, 



un = 0, n • [Vu + (Vu) T ] • (I - fin) = 0, (r g Sf), 



(2.3) 



(2.4) 



where I is the metric tensor. 

In the porous medium, the flow is driven by the capillary pressure in the menisci that collectively form 
the wetting front, and, for the problem in question, the characteristic pressure and velocity are 2a /a and U, 
respectively. Using the notation p and u for the dimensionless pressure and velocity in O2 , in a frame moving 
with the contact line we have an equation of motion in the form of Darcy's law, 



u-U = -iCVp, (ren 2 ), 



(2.5) 



where K — 2a k / '(fiaLU) is the nondimensionalized permittivity of the porous matrix (k is the actual permit- 
tivity) and U = XJ/U is a unit vector directed along the velocity of the porous substrate. Given that in porous 



media k oc a 2 (Probstein 



1989), we have that K = 0(e/Ca) and, for the problem to be nontrivial, it is assumed 
to be finite in the limits (|2. If) . f|2 . 2|) . In the present context, it is convenient to define L by setting K = 1. 
Since the porosity <fi in a homogeneous matrix is constant, the mass balance equation has the standard form 



V • u = 0, so that, after substituting (|2.5p into it, one arrives at Laplace's equation for p in f2 2 : 



V 2 p = 0, (re n 3 ) 



(2.6) 



On the wetting front, one has the kinematic condition that the front propagates with the velocity of the fluid, 
i.e. in the reference frame described above where the process is steady, 



u • n = 0, (r 6 S w ), 



(2.7) 
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and needs to specify the dynamic condition on the pressure p. To formulate this condition in a general case, 

one has to consid er the modes of motion the men isci go through on the pore scale as the macroscopic wetting 



front propagates (S hikhmurzaev fc Sprittles 



2012t ?). Here we will be considering the simplest case involving 
the main, 'wetting', mode. In this case, for an 'effective' pore with a circular cross-section the dimensionless 
pressure at the front is given by 

p = -cos(9 d , (2.8) 
where 9 d is the contact angle formed on the pore scale by the repr esentative meniscus with the pore wall 



(Fig. [T]) . Unlike the simplified approach pioneered by 



Washburn! (|192l[ ). where it is assumed that 6 d is constant 



and equal to the prescribed static co ntact angle 9 S , we ne ed to take into account that, as demonstrated by 



Shikhmurzaev 



20071 for a review), 9d depends on the wetting speed, i.e. 



numerous experiment (e.g. see Ch. 3 of 
to consider the dependence 

(u-V)-n = Uy(6 d ), (reS w ), (2.9) 

where XJ* cl is the appropriate velocity scale depending on the material parameters of the system (as all velocities 
above, we have it nondimensionalized using U). The function f(9 d ) has to be determined theoretically or 
empirically. For ex ample, one can appl y f(9 d ) derived using the theory of dynamic wetting as a process of 



interface formation ( Shikhmurzaev 



2007), which has been shown to reliably describe experimental data, though, 



in the present context, any function f(9 d ) representing the experimentally observed dependence of the form 



(Blake & Havnes 


— , „ 
1969 


). For more 


recent review 1 


Shikhmurzaev 


2011) 



1969 ). For more information about the dynamic wetting modelling we refer the reader to a 



From the theory of flows with forming interfaces one has 

m v = / (! + (!- pj e )cosfl 8 )(cosfl B - cosfl d ) 2 x 1 2 
S( d) \ A(cos9 s + B)(cos9 d + B) 

where B = (1 - pfj-^l + p s le u (9 d )), 



m _ u cl _ i (iPi 0) (i + ^py 

U '] = 



1/2 



u u 



(2.10) 



u (9 d ) 



sin V d - v d cos ti d 



sin ti d cos v d - v d 

and pf e , a, /?, 7, r are m aterial constants characterizing the contacting m edia. Their values for some 



systems can be found elsewhere ( Shikhmurzaev 



2007 



Blake k, Shikhmurzaev 



2002). 



So far, the flow in the pure fluid and in the porous medium were considered separately, and to link them one 
has to specify three boundary conditions at Sq. One condition that we obviously have on this surface is the 
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continuity of mass flux: 

(u - <£u) • n = 0, (re So)- (2.11) 

For the velocity components parallel to So, i.e. for ii|| = u ■ ( I — no no ) and U|| = u ■ (I — npno), a number of 



boundary condit i ons have been discussed in the literature (e.s. ISaffman 



1999; 



Nield 



2009; 



Auriault 



1971 



Jones 



1973 



20101) , following the experiments repo rted by 



Auriault 



20ll 



Mure och & Soliman 



Beavers fc Joseph! (jl967h and an empir 



2010), these conditions are aiming at 



ical condition these authors proposed. As noted, for example, by 
capturing the effects of order 0(e), i.e. go beyond the classical, i.e. Oth-order, approximation in the continuum 
limit (|2.1I) . In the 0th approximation, all these conditions reduce to no-slip for the pure fluid, 



U||=U||, (re S ), 

and it is this condition that we will be using here. 

The condition of continuity of pressure on So in the dimensionless form yields that 



(2.12) 



P = P- 



2 ' 



and hence, to leading order in the limits (|2.ip and (|2.2I) . one has 



p = 0, (re So). 



(2.13) 



Importantly, unlike the case of an impermeable solid substrate, where one has to specify the dynamic contact 
angle (as we need to specify 8d on the pore scale), here, for the effective contact angle 6q, we will require only 
that the flow parameters in the porous medium remain regular at the contact line and that the solution in the 
pure fluid exists. 



3. Dynamic contact angle 

Consider the asymptotic behaviour of the solution to (|2.3[) - (|2.13p in the case of a two-dimensional flow as 
the distance to the contact line r — > 0. In the polar coordinates (r, 9) shown in Fig. [TJ to leading order, for the 
pressure p in the porous medium one has 

d 2 P i d P i d 2 P / a / m 

or 1 r or r z o0 z 
dp 

p(r,0) = 0, —(r,-6 w ) = rsm0 w , 
where the last condition follows from (|2.7p and (|2.5I) . The separable solution to this problem is obviously given 

by 

p = rtan^u, sin#, (3-1) 
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so that, using the dynamic boundary condition (12.81) . we obtain that 

Then, from (|2.9I) , where now (u — U) • n = sin 9 W , we have an equation determining 8 W : 

sm9 w = U^f(7r/2). (3.2) 
In order to consider the flow in the pure fluid, we introduce a stream function 



1 dip _ dip 
r oO or 



(3.3) 



so that for the leading-order term ijji of the asymptotic expansion of ip as r — > we have a biharmonic equation 

T^ + -V + -^w ^= ' (O<r,O<0<0D), (3.4) 
or 2 r or H d8 z I 



together with conditions (|2.4p . i.e 



90 

on the free surface, and conditions (|2.1ip and (|2.I2D . i.e. 



Mr,0 D ) = 0, -Jt±( r ,9 D ) = 0, (3.5) 



^(r,0) = r^tan^, ^( r ,0) = r, (3.6) 

on the surface of the solid substrate. In writing down the first of conditions Q3.6P we made use of (|3.1[) and 
integrated along So- 

The separable solution to the problem (|3.4p ~ (|3.6p has the form 

ipi = r(A x sin6> + A 2 0sin0 + A 3 cosd + A 4 dcosd), (3.7) 



where the constants A\, . . . ,A± are given by 



e D cos 2 e D 

A l = — — — n ~7i n A 3~ 



sin 9d cos 6d — Od sin^^ cos6*£> — 9d 



A, = sin °d + A sin6> D cos 8 D 

sin^D cos(?d — 9d sin#D cos#d — 8d ' 

A 3 = 0tan6»„,, A 4 = A 2 cot 6>_c . (3.9) 

Using (|3.3p . we can write down the radial projection of the second equation in (|2.3D in the form 

dp _ n d 3 + 1 d 3 + 1 d 2 \^ 
dr \r dr 2 d6 r 3 d6 3 r 2 drdB ) ' 

and, after substituting the solution (|3.7I) . arrive at 

dp 2 

— = — ^(sin# + cot^D cos 9)A 2 - 
dr r A 



Thus, the leading term in the coordinate expansion of the stream function will not give rise to a nonintegrable 
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Figure 2. Dependence of the effective dynamic contact angle 9c on the dimensionless contact-line speed U» for 
different porosities. Curves 1, 2, 3 and 4 correspond to cj> = 0.1, 0.2, 0.3 and 0.4, respectively. 



singularity of pressure if and only if A 2 — 0, or, using (|3.8I) and (13.91) . 



tan#n = — </>tan( 



Given the expression p.2[) for sin6 w and introducing the dimensionless contact-line speed E/» = l/(/(7r/2)C/* z ) 
U / {f (it / 2)U ci) : we can also write this equation down as 



Id = 7r + arctan — - 



(3.10) 



This equation specifies Op in terms of the speed U of the contact line with respect to the substrate, the porosity 
of the matrix <j> and the material properties of the contacting media accumulated in U c i and f(jr/2). The 
velocity-dependence of 6n for different porosities is shown in Fig. [3] As one can see, this dependence is much 
steeper than the velocity-d ependence of the dyn amic contact angle for the dynamic wetting of an impermeable 



substrate (e.g. see Ch. 3 of 



Shikhmurzaev 



2007, for a review). In particular, d6r>/dU* — > +oo as [/* — s- 1+. 
The dynamic contact angle given by (|3.10l) ensures integrability of the normal stress on the free surface and 
hence the existence of a solution. The next term in the asymptotic expansion of ip has the form ip2 = r 2 F(6) 
and therefore can give rise, at most, to a logarithmic (i.e. integrable) singularity of p at the contact line, which 
does not affect the existence of the solution. 

The obtained result has a clear physical meaning. In the case we are considering, the dynamics of imbibition 
determines the normal to the substrate component of the pure fluid's velocity independently of the pure fluid's 
bulk flow, and, given that the tangential component of the fluid's velocity on the surface of the substrate is 
prescribed, as it satisfies the no-slip condition with a known speed of the substrate, we have a moving contact- 
line problem where on the solid surface both components of velocity are set. In this situation, in the pure 
fluid the solution exists only if from the family of stream functions described by (13.7[) we choose the one that 
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corresponds to a uniform flow (A2 = 0, A4 = 0) with the projections of velocity on the normal and tangential 

to the substrate directions equal to the speed of imbibition and the speed of the solid substrate, respectively. 

The solution (|3.10[) for a steady spreading of a fluid over an unsaturated substrate ceases to exist when the 

dimensionless velocity U* becomes equal to 1 and both 9r> and W reach tt/2. For C7* < 1 we have a different 

regime with the two contact lines CL1 and CL2 no longer coinciding as the wetting front moves ahead (dashed 

line in Fig.[I|, and the pure fluid finds itself on a saturated substrate. Then, in the vicinity of CL1 the pressure 

p satisfies Laplace's equation V 2 p = together with the boundary conditions 

p(r,0) = 0, g(r,-7r) = 0, 

and the local solution with the regular pressure gradient has the form 

P = fjC„r 1 /2+« sin Q +n ^ ; (3.11) 

where C n , n = 1,2,... are constants. Then, according to (|2.5I) and (|2.1ip . one has ug(r, 0) oc r 1 / 2 — >• as r — > 0. 
As a result, if u r (r, 0) = 0(1) as r — > 0, then, to leading order as r — > 0, we will have the 'classical' moving 
contact-line problem, with no i mbibition and the n o-slip boundary condition on the solid. This problem, as is 
well known (see, e.g., Ch. 3 of 



Shikhmurzaev 



2007), has no solution. The only way out of this situation is to 
conclude that, as CL1 and CL2 separate, the contact line CL1 stops moving. Then, for the pure fluid one simply 
has a static contact line with the imbibition velocity oc r 1 / 2 near it and Qd = tt/2. It is easy to verify that the 
solution to this problem exists. As CL1 and CL2 separate and CL2 moves ahead, the imbibition process near 
CL1 slows down, leading to C\ = in p. lip and hence ug oc r 3 / 2 , so that 6d no longer has to be equal to 7r/2 
and can go down as the imbibition continues. 



4. Discussion 



The described scenario is in agreement with the available experimental data. In experiments on the spreading 
drops, it has been observed that the regime where the drop base and the saturated area underneath it expand 
together is followed by the regime where the two conta ct lines, CL1 and CL2, separate, as drop's base stops 



expanding whereas the satur ated area continues to grow (jClarke et al 



2007 



Markicevic et al 



2002 



Starov et al 



2003 



Keshav fc Basu 



2010J). This is usually attributed to the 'competition' between spreading and imbibition, 



a nd the above analysis sh ows what this actually means. 



Markicevic et al 



(2010) report that in their experiment the spreading droplet maintained the shape of a 



spherical cap throughout the process and that its base stopped expanding when "the droplet is a half of 
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sphere". In other words, it stopped expanding when the dynamic contact angle 9d became equal to 90°. This 
is exactly what follows from the result of our analysis. 

The dependence of the effective contact angle 9d on the contact-line speed for U* > 1 given by (|3.10p and 
illustrated in Fig. [2] awaits its experimental verification. The issue here is that, to extract this dependence 
from experiments w i th un steady flows, one has to deal with very short time intervals (<C 40 ms according 



to 



Markiccvic et al 



( 20101 )) with a well-controlled spatial resolution and conditions of the Darcy-scale de- 
scription satisfied. The problem is made more complicated by the fact that, as one can see in Fig. [2 the 
velocity-dependence of 9d is very steep, which brings in additional conditions on the temporal resolution of the 
experiments. On the other hand, experiments with steady flows of the kind commonly used to study dynamic 
wetting of impermeable solid surface require sizeable substrates with well-reproducible properties. To date, no 
systematic measurements of 9d with controlled spatial resolution have been reported. 

As we have shown, the requirement of the absence of a nonintegrable singularity of pressure in the pure fluid 
at the contact line, which is equivalent to the requirement of the existence of a solution, uniquely specifies the 
velocity-dependence of the 'effective' dynamic contact angle 9d formed by the free surface and the 'effective' 
surface of the porous substrate. In practical computations, it may be convenient to set (|3.10|) as a boundary 
condition while considering the normal-stress condition on the free surface as the equation determining the free 
surface shape. One implication of the obtained result is that, if one specifies 9d as an additional boundary 
condition different from (|3.10p . then, for a solution to exist, it will become necessary to use a slip boundary 



2000). As mentioned earlier, 



condition instead of the no-slip condition (I2.12[) employed here ([Davis fc Hockind l: 
this would mean bringing in effects of 0(e) into the essentially classical, i.e. O(l) as e — > 0, formulation. If 
the velocity-de pendence of 9p different fr om (13.101) is impose d together with the no-slip condition (|2.12[) for 



the pure fluid 



Alleborn &: Raszillier 



2004 



Reis Jr. et al 



2004 ). then there will be no solution to the problem, 



although this fact can be hidden behind, and as a result masked by, simplifying assumptions made in the process 
of finding the solution (e.g. lubrication approximation) and the numerical implementation. 

From the theoretical viewpoint, it would be interesting to consider the essentially unsteady process of sepa- 
ration of CL1 and CL2, which marks the transition between the two regimes described above. As with every 
finite-time transition, this is a challenging problem deserving a detailed investigation. 
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